Stored Electromagnetic Energy and Antenna Q 

Mats Gustafsson*and B. L. G. Jonsson 1 ' 
November 26, 2012 



Abstract 

Decomposition of the electromagnetic energy into its stored and radiated parts 
is instrumental in the evaluation of antenna Q and the corresponding fundamental 
limitations on antennas. This decomposition is not unique and there several proposals 
in the literature. Here, it is shown that stored energy defined from the difference 
between the energy density and the far field energy equals the new energy expressions 
proposed by Vandenbosch for many cases. This also explains the observed cases with 
negative stored energy and suggests a possible remedy to them. The results are 
compared with the classical explicit expressions for spherical regions. It is shown that 
the results only differ by a factor ka that is interpreted as the far-field energy in the 
interior of the sphere. Numerical results of the Q-factors for dipole, loop, and inverted 
L-antennas are also compared with estimates from circuit models and differentiation 
of the impedance. 

1 Introduction 

It is well known that the electrostatic energy in free space can be written as an integral 
of the energy density, eo\E\ 2 /A, or equivalently as an integral of the electric potential, <fi, 
times the charge density, p, [1-4]. A similar expression holds for the magnetostatic energy. 
The electrodynamic case is more involved. In [5], Carpenter suggests a generalization in 
the time domain based on (f>p + A- J, i.e., the sum of the scalar potential times the charge 
density and the vector potential, A, times the electric current density, J, see also [6,7]. 
In [8] , Vandenbosch presents general integral expressions in the electric current density for 
the stored electric and magnetic energies. These expressions are similar to the expressions 
by Carpenter but include some correction terms, see also [9]. 

The expressions by Vandenbosch are very useful to analyze small antennas [10-13] 
and have been verified for wire antennas in [14]. One minor problem with the proposed 
expressions is that they can produce negative values of stored energy for electrically large 
structures [11]. On the other hand, the classical work by Chu [15] is based on subtraction of 
the power flow and explicit calculations using mode expansions of the stored energy outside 
a sphere, see also [16]. This gives simple expressions for the minimum Q of small spherical 
antennas [15,16]. The major shortcoming is that the results are restricted to spherical 
regions although some generalizations to spheroidal regions are suggested in [17, 18]. The 
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results have also been generalized to the case with electric current sheets by Thai [19] and 
Hansen and Collin [20] by adding the stored energy in the interior of the sphere. 

In this paper, we investigate stored electric and magnetic energy expressions based on 
subtraction of the far-field energy density. The expressions are suitable for antenna Q and 
bandwidth calculations and closely related to the classical methods in [16] and others for 
antenna Q calculations. They are not restricted to spherical geometries and, furthermore, 
resembles the recently proposed expressions by Vandenbosch [8]. The results provide a 
new interpretation of Vandenbosch's expressions [8] and explain the observed cases with 
negative stored energy [11]. They also suggest a possible remedy to the negative energy 
and that the computed Q has an uncertainty of the order ka, where a is the radius of the 
smallest circumscribing sphere and k the wavenumber. This is consistent with the use of 
Q for small (sub wavelength) antennas, where ka is small and Q is large [15, 16]. Analytic 
results for spherical structures show that the expressions in [8] for Q differ with ka from 
the results in [20], that is here interpreted as the far-field energy in the interior of the 
sphere. The results for Q are also compared with estimated values from circuit models 
and differentiation of the impedance [21, 22] for dipole, loop, and inverted L antennas. 

The paper is organized as follows. In Sec. 2, the stored electric and magnetic energies 
defined be subtraction of far-field from the energy density are analyzed. Analytic results 
for spherical geometries and comparison with classical results are presented in Sec. 3. The 
coordinate dependence is analyzed in Sec. 4. Stored energies from small structures are 
derived in Sec. 5. Comparisons with numerical results for dipole, loop, and inverted L 
antennas are given in Sec. 6. The paper is concluded in Sec. 7. 



2 Stored electromagnetic energy 

We consider time-harmonic electric and magnetic fields, E(r) and H(r), respectively, 
with a suppressed e~ lw * dependence, where co denotes the angular frequency. The Maxwell 
equations in free space are [1] 

V x E = iuj/j, H = irj kH 

V x H = -ioje E + J = -^E + J, 

where J denotes the current density, while eo, po, and 770 = y^o Ao are the free space 
permittivity, permeability, and impedance, respectively. For simplicity, we interchange 
between the angular frequency and the free space wavenumber k = w/co, where the speed 
of light Co = l/^/^o e o- We also recall the continuity equation, V • J = iuip, relating the 
current density J with the charge density p. 

It is widely accepted [1, 2, 4] that the time-harmonic electric and magnetic energy den- 
sities are eo\E\ 2 /A and pq\H\ 2 /4:, respectively. On the other hand, there are a few alterna- 
tive suggestions in the literature [3], and the energy densities are not observable [5]. The 
electric and magnetic energies comprise both radiated and stored energies; however, for 
antenna Q calculations one must distill the stored energy. In this section we analyze stored 
electric and magnetic energy expressions suitable for antenna Q and bandwidth calcula- 
tions, which shed new light on the meaning and practical applicability of the methodology 
for evaluating stored energies directly from the sources that was derived recently in [8] . In 
subsequent sections of this paper we elaborate the implications and practical applicability 
of the results of this section. 
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Figure 1: Illustration of the object geometry V and with outward normal unit vector n 
and current density J(r). The object is circumscribed by a sphere with radius a. 



It follows from Maxwell's equations (1) that the sources and fields obey the conservation 
of energy equation in differential form, 

i2 W (||£| 2 - ^\Hf) + 1 -E.r = ^V.(Ex H% (2) 

where the superscript * denotes complex conjugate. We consider current distributions J 
whose support is in a volume V bounded by the surface dV, see Fig. 1. Integrating (2) 
over this volume gives the real part result 

^ / E(r)xH*(r).h(r)dS = ~ [ J3(r) • J*(r) dV, (3) 
1 Jav 1 Jv 

where h denotes the outward- normal unit vector of the surface dV. The first term in the 
real part expression (3) is readily identified in view of the Poynting vector [1,4] as the 
time-average radiated power flow through the surface dV, so that (3) equates the radiated 
power exiting dV to the time average of the power generated by J, as expected from 
energy conservation. Furthermore, integrating (2) over all space shows that the radiated 
power exiting the surface dV can be expressed in terms of the far field as 

Pi = Re / ^)xjr W . ftWds= nmp?^ , 

Jav 2 Jn 2 Vo 

where O denotes the surface of the unit sphere and the far field behaves like E(r) ~ 
e lkr F(r)/r as r — > oo, where r = rr and r = |r|. Similarly, by integrating (2) over all 
space one obtains the imaginary part result 

£ > W f-^>MfdV = Im / E(r) ; J '"-> dV, (5) 



4 4 Jy A.UJ 

where we used the fact that the integral of the imaginary part of the divergence term in (2) 
vanishes as the integration volume approaches M 3 . The imaginary part result (5) relates 
the well-defined difference between the time-average electric and magnetic energies with 
the net reactive power delivered by J. 

As is well known [15, 16], the total energy, defined as the integral of the energy density 
integrated over all space, is unbounded due to the 1/r 2 decay of the energy density in the 
far radiation zone. This is resolved by decomposition of the total energy into radiated and 
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stored energy. The stored energy is, however, difficult to define and interpret. The classical 
approach used by Chu [15] and Collin & Rothschild [16], and subsequently by others, is 
based on mode expansions, and therefore restricted to canonical geometries. Spherical 
regions are most commonly considered but there are also some results for cylindrical [16] 
and spheroidal [17, 18] structures. The stored energy density is customarily defined as 
the difference between the total energy density and the radiated power flow in the radial 
direction, thus the stored electric energy becomes 

wP = ^f \E(r)\ 2 -r, Re{E(r)xH*(r)}-rdV, (6) 

where the subscript r in Mj? = {r : ]im r0 -^. oo |r| < ro} is used to indicate that the integration 
is over an infinite spherical volume. The classical results by Chu [15] are for spheres with 
vanishing interior field [16], so that the stored energy is due to the exterior field only (i.e., 
for the region where r > a where a is the radius of the smallest sphere circumscribing 
the sources). The Thai bound [19] generalizes the results to fields generated by electric 
surface currents, see also [20]. Here it is observed that there is a stored energy but 
no radiated energy flux in the interior of the sphere. The definition (6) is useful for 
spherical geometries and can be generalized to cylindrical geometries [16] . It is difficult to 
generalize it to arbitrary geometries due to its coordinate dependence that originates from 
the scalar multiplication with r. The subtraction of the radiated energy flow is equivalent 
to subtraction of the energy of the far field outside a circumscribing sphere, c/., (4). This 
suggests an alternative stored electric energy defined by subtraction of the far-field energy, 
t. e. 

<' = f/ i; | £W P-E^dV, (7) 

where the integration is over the infinite sphere Rj?- 

We note that the definitions with the power flow (6) and far field (7) differ only in the 
interior of the smallest circumscribing sphere associated with the source support. In the 
interior of the smallest circumscribing sphere, which we assume next to be of radius a, 
this subtracted far-field energy is then 



eo 
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£ JjF(r)\ 2 <mdr = £-P T . (8) 



Assuming that the contribution to the true stored electric energy, say We, due to the 

(p) 

exterior field outside the smallest circumscribing sphere, is equal to that of W E and 

wP in (6) and (7), and that it subtracts some non-negative value less than eo\F\ 2 /(Ar 2 ) 
inside the sphere, then we obtain the bound 

WP <W E < WP + ^-P T . (9) 

ZCo 

This means that the stored electric energy can be bounded from below and above by (7). 
The stored magnetic energy, Wyp , is defined analogously. It is customary to normalize 
the stored energy with the radiated power to define Q- factors. The Q- factor is Q = 
ui8cx{Qe,Qm}, where 

Qe = —5 — and Q M = — 5 — (10) 
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and we have included a factor of 2 in the definitions of Qe and Qm to simplify the 
comparison with antenna Q. This translates the bound (9) into 

max{0, Q (F) } < Q < Q (F) + ka, (11) 

where we have added that Q is non-negative. 

We show that the stored energy with the subtracted far field (6) is similar to the energy 
defined by Vandenbosch in [8] for the vacuum case. For simplicity we express the energy 
using the scalar potential <p and the vector potential A in the Lorentz gauge [1,2,4], so 
that (V 2 + k 2 )(f>(r) = -p(r)/e and (V 2 + k 2 )A(r) = -fi J(r) and therefore 

4 > {r) = e \G*p){r) = - I G(\r - r'\)p{r') dV' (12) 

e o Jv 

and 

A{r) = p (G * J)(r) = p [ G{\r - r'\)J(r') dV', (13) 

Jv 

where * denotes convolution and G is the outgoing Green's function i.e., G{r) = e lkr / (4-7rr). 
The vector and scalar potentials are related by V-A = ifc^/co and the electric and magnetic 
fields are given by [1] 

E = \uA - V<j) and H = p^ 1 V x A. (14) 
We also use the corresponding far-field potentials defined by 

<Poo(r) = ^ r J v P(r')e' ik ^'dV' (15) 

and 

A 00 (r) = f^- f J(r')e- ifefV dV' (16) 

47T Jv 

giving the electric far-field 

F(r) = iLuA 00 (r)-rik(f) 00 {r). (17) 

Using that the far-field is orthogonal to f, i.e., r ■ F = 0, the far-field radiation pattern 
obeys 

l^^^iAoo^)! 2 -^!^^)! 2 . (is) 

The electric energy density is proportional to 

\E\ 2 = w 2 | A\ 2 - 2 Re{iwA • V<f } + | V0| 2 

= u 2 \A\ 2 - 2k 2 \0\ 2 + |V0| 2 - 2Re{iwV • (<f A)}, (19) 

where we used V • (4>* A) = cp*V ■ A + A ■ V0* = \k\4>\ 2 / c + A ■ V(p*. We integrate this 
result over a large sphere to get the far- field type stored electric energy (7) expressed in 
the potentials 



4W£> 



f \ E (r)\ 2 - l -^^dV= [ \V<P(r)\ 2 - k 2 \Hr)f 



+ u 2 [ \A(r)\ 2 - \ A °°(pl) _ k 2 ( \ Hr) \2 _ dV , (20) 
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where we applied the divergence theorem to the integration of the last term in (19), 
obtaining via the discussion in (17) and (18) that 

/ lm{(p* (rr)A,.(rr)}r 2 dn -»• (21) 
Jn 

as r — > oo, see (17). 

Use the energy identity for the Helmholtz equation, | V0| 2 — /c 2 |^| 2 = eg 1 Re{0p*} + V • 
(Re{</>*V^>}), and that (jfVcj) — > \kr\(j)\ 2 for large enough r, to rewrite the first two terms 
in (20) as 



/ 

Jm 



\V(p(r)\ 2 -k 2 \<p(r)\ 2 dV = e 1 Re / <f>(r) p* (r) dV 



v 



= / / rtrO ^-pjV MdV.dV,, (22) 
JvJv 4irei\ri-r 2 \ 



where we also used that the surface term vanishes. The Green's function identity, see 
App. A 



/ 

it 



G { \r - r ll)0 ' (i r - „|) - dV _ + ,23) 



where ji(.z) = (sin(z) — zcos(z))/z 2 is a spherical Bessel function [4], is used to rewrite 
the two remaining terms in (20) as 



r r 2 2 

+ i/ / J(r 1 )-%^ 2 j 1 (fc|r 1 -r 2 |)J*(r 2 )dV 1 dV 2 (24) 
JvJv ovrri2 



and 

-ikr'-r n ( m l\ A\T r |2 

dV = - 

f f . . r- 

+ i 



/ 



2 \f v e- [kr ^(OdVT „. /■/■ , . sin^ln-^l) 
|G * p| - 16 ,2 r 2 dV = - / / P(ri ) ^ p (r 2 ) dV, dV 2 



CP 2 2 

1 / P^O^^ji^ln-rsD^Md^d^. (25) 



We note that the first terms in the right-hand side of (24) and (25) only depend on the 
distance n 2 = — **2| and are hence coordinate independent, whereas the last terms 
depend on the coordinate system due to the factor r\ — r\ = (rq + r 2 ) • (rq — r 2 ). The 
coordinate dependence originates in the explicit evaluation of the integral in (23) over 
large spherical volumes Mj? that is necessary due to the slow convergence of the integral 
in (23), see also App. A. 

Collecting the terms in (22), (24), and (25), we get a quadratic form in the current 
density J for the far-field type stored electric energy (20) as 

Wf = <°) + <J) + < M *), (26) 
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where wj^ ^ + is the coordinate independent part 

<"> + *£}> - ? JJ v v. ■ JW, ■ rM^f 

- {k 2 J(n) • J(r 2 )* - Vi • J(n)V 2 • J*(r 2 )) Sm(fc| ^~ r2|) dV x dV 2 (27) 

and VF"g F °' ) and W^m^ contains the cos and sin parts, respectively. The coordinate depen- 
dent part is 

^EM 2) = f JJ v Im {k 2 Ji ■ J*2 ~ Vi • JiV 2 • ^2} gj^^i Ji(fcki ~ r 2 |) dV, dV 2 , (28) 

where J n = J(r n ), n = 1,2. The coordinate independent part Wg' ' + W^m^ is identical 
to the energy by Vandenbosch in [8] for vacuum and hence presents a clear interpretation 
of the energy [8] in terms of (7). We also see that the definition (7) explains the peculiar 
effects of negative stored energies [11] and suggests a remedy to it in (9). The coordinate 
dependent part W^ff is more involved. Obviously the actual stored energy, as any physical 

quantity, should be independent of the coordinate system. First, we observe that W^ff = 
for any current density that has a constant phase. This includes the fields originating 
from single spherical modes on spherical surfaces and hence most cases in [15, 16, 19,20]. 
It also includes currents in the form of single characteristic modes [12]. We also get the 
coordinate independent part by taking the average of the stored energy from J and J*. 
The term wj^Ji is further analyzed in Sees 4 and 5. 

For the stored magnetic energy we can use |S| 2 = |V x A\ 2 or simpler the energy 
identity (5), to directly get the difference 



/ 

Jm 



Li \H\ 2 - e \E\ 2 dV = Re [ A ■ J* - <f>p* dV, (29) 

Jv 



IV 

where we used 

E ■ J* = koA ■ J* - V ■ (</>J*) - ]u<f>p*. (30) 

This gives the far-field type stored magnetic energy wffi = + W^j' + ^em^ 

where the coordinate independent part 

(Fo) , w (Fi) _ Mo f f 7/ \ T */ .cos(fc|n -r 2 |) 



< 0) +wfi ) = ^/ y / v J(n)-J*(r 2 ) 



47r|ri — r 2 | 



- (k 2 J(n) • J(r 2 )* - Vi • J(n)V 2 • J*(r 2 )) r2|) dV, dV 2 

is expressed as a quadratic form in J, see also [8]. We also have the radiated power 



(31) 



Pr = I Ivlv ( ^ J(n) " J(ra) * " Vl " J(n)V2 " r{r2)) ^\ri- rf dVl ^ • (32) 

It is illustrative to rewrite the coordinate independent far-field stored energy in the 
potentials: 

r (F ) _ Re f „ , „,fF„i Re 



^ " "4 



jTpVdV, <f o) = ^^J*-AdV (33) 
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and 



w (Fi) Re [ k f T* dA *^jv CiA\ 

where it is assumed that the frequency derivative of J and p are negligible in (34). We 
note that the sum of the first terms, W^ ^ + wj^°\ corresponds to a frequency- domain 
version of the energy expression by Carpenter [5], see also [6, 7]. Moreover, they reduce to 
well-known electrostatic and magnetostatic expressions in the low-frequency limit [1] . 

It is also convenient to follow standard notation in the method of moments (MoM) 
and introduce the operators £ e and Cm such that C = C c — Cm is the integral operator 
associated with the electric field integral equation (EFIE) [23]. Here, the operators are 
generalized to volumes and defined from 

(J, £cJ} = ^ J J Vi • J(n)V 2 • J*(r 2 )G(|n - r 2 \)dY 1 dV 2 , (35) 

{J,CmJ)=ik[ I J(n)- J*(r 2 )G(|ri-r 2 |)dV 1 dV 2 , (36) 
JvJv 



and 



{J,CemJ) = %l I CjVi- J(ri)V 2 - J*(r 2 ) 



2 .lv Jv ^ 



-»(„)■ J. W ) gg^p!) dV,dV,. (37) 



They are defined such that the stored electric and magnetic energies and radiated power 
are 

Wf 0) = ^Im(J,£ c J) (38) 
^ o) = £lm(J,£ m J> (39) 
W^^g-Im^^J) (40) 

Pr= T jRc(J,(C c -Cm)J}- (41) 

Efficient evaluation of the C operator is instrumental in MoM implementations where 
the discretized versions are often referred to as impedance matrices. The relations above 
show that the corresponding matrices for the coordinate independent stored and radiated 
energies are available by evaluating the real and imaginary parts of the MoM impedance 
matrices with the addition of the mixed part (37). The stored energy is hence computed 
with negligible additional computational cost in MoM implementations. Moreover, (41) 
shows that Re C is positive semidefinite. 



3 Electric surface currents on a sphere 

The two formulations (6) and (7) for the stored energy can be compared for electric sur- 
face currents on spherical shells. This is the case analyzed by Thai [19] and Hansen & 
Collin [20]. We expand the surface current on a sphere with radius a in vector spher- 
ical harmonics Y, see App. B. For simplicity, consider the surface current J(r) = 
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•A) Y T <Tml(r)S(r — a). It induces the electric and magnetic fields 

E(r) = ivoJo / and ff(r) = J ^ / , (42) 

(p) (v) 
where p = 1 for r < a and p = 3 for r > a, 2 is the spherical vector waves, and 

the radial functions in Hansen [24], defined as 

>\k) r = l 

2 



R#(k) = < I'a^w) ( 43 ) 



5k 

where z, 1 ^ = are Bessel functions, = n^ Neumann functions, = hf 1 ^ Hankel 

functions [24], and k = ka. We note that the derivatives of R^(k) are easily expressed in 
z^ p \ see App. B. Here, r = 1 is transverse electric (TE) and r = 2 transverse magnetic 
(TM) waves. Moreover, the dual index fisr = 2ifr = l and f = 1 if r = 2. The 
current in (42) is rescaled as Jo = Jo R^i{ka) (ka) and below we let Jo be real valued 
to simplify the notation. We also note that the coordinate dependent term (28) vanishes 
for single spherical modes. 

3.1 Far-field type stored energy for the TE case 

We start with the transverse electric (TE) case r = 1, i.e., J(r) = Yiami(r)5(r — a) 
that is divergence free, V • J = 0. The integrals in (26) are evaluated analytical by 
expanding the Green's functions in (35), (36), and (37) in spherical modes, see App. B. 
Using V • Yicrmi = 0, we get (J, £ c J) = for (35) and hence the first part of the stored 
electric energy H^ F °' = 0. The expansion of the full Green's dyadic, G = GI, (86) gives 

( J > £m J ) = j v j v Yiomj(ri)*(ri - a) ■ G(|n - r 2 |) ■ Y l(Tm i(r 2 )5(r 2 - a) dV x dV 2 



= a 4 / / Yi CTm Kr-i)-G(|r 1 -r 2 |)-Yi CTm Kr 2 )d^id^ 2 = ia 4 A ; R^( K )R^(K) (44) 
JnJn 

for the terms in (36) to get the first part of the stored magnetic energy from (39) 
as 4wr/ " 1 W I [J ?o) = -aW^R^Rjf. The radiated power follow from (41) 2r/ " 1 P r = 

— Re( J, Cm J) = cl 2 k 2 Jq (R^) 2 . The corresponding expansion of the frequency derivative 
of the Green's function (86) is used for the terms related to (37) 

77^72 (J,C cm J)= [ [ Yw(f i) • d Q(|ri ~ r2|) • YWf 2) dOi dQ 2 
\k A a A J£ J n J n dk 

= i iL («rS ) (' s )rS ) («)) =K«rS ) («)4 1) («))' 

= i(RS?RSl ) + .RSf'RSl ) + .RS ) R«0, (45) 

where ' denotes differentiation with respect to k, giving 4lo7]q 1 W^^ = —^^-Jq(k R^ R^)'- 
Collecting the terms gives the electric and magnetic Q-factors as 

(F) _ 2u,wf >(«) _ (^R^RgW f4fi x 

Qi '- e( } " ~^r" - 2(r«(«) )2 (46) 
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and 

RiJ ; (k) 

respectively. We note that = j; and = n^ can be used to rewrite the Q-factors, 
however the form with the radial functions simplifies the comparison with the TM case 
below. The differentiated terms are also easily evaluated using (45) and (84). 

3.2 Far-field type stored energy for the TM case 

The transverse magnetic (TM) case is given by r = 2 and generated by the current density 
J(r) = Jo Y2ami(r)S(r-a) that has the divergence V-Y^mi = + 1) Yami /r. With 

the expansion of the Green's function (85) we get the part related to the charge density (35) 

^(J,£cJ)= / / V 1 -Y2and(r 1 )G(\r 1 -r 2 \)V2-Y2and(r 2 )dn 1 dn 2 
a Jq JnJn 

^iWC'W (48) 

and the full Green's Dyadic expansion (86) gives 

7^2<J,£mJ>= f j Y2ami(ri)-G(|ri-r 2 |)-Y:w(f2)dnidn2 

= I* (bS)(«)bSPW+ 1(1 + 1)^^V (49) 

for the part related to the current density (36). The expansions of the frequency derivatives 
of the Green's function (85) and Green's Dyadic (86) give 

- Vi • Y 2CTm Kn) 8G(l fc r 2 g~ r2|) V 2 • Y2omj(r 2 ) dfii dft 2 

= 21(1 + l)n t (K) hi*) ~ k 2 (krW(k) R 2 f{K))'. (50) 

for the part related to (37). 

Collecting the terms gives that the normalized radiated power is 

2 ??0 - 1 P r / J 2 = Re(J, (Co - C m )J)/Jo = a 3 K(R ( V) 2 . (51) 
The electric and magnetic Q factors are finally determined to 

(F,,, M'i'MRg'M)' 

Q "* W - 2(Rg' W ). <52) 

and ^ 

QL,M = Q^eM (i) ' ( 53 ) 

R 2/ (k) 

respectively. We note that the expressions for the TE case in (46) and (47) and TM case 
in (52) and (53) are written in identical forms by using the radial functions (43). 
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3.3 Power flow stored energy W^' 



EM 



The stored electric energy with the subtracted power flow (6) is analyzed by Hansen & 
Collin [20], see also Thai [19]. The integral (6) is decomposed into integration of the exte- 
rior and interior regions where we have outgoing waves, ,, and regular waves, u~mZ' 
respectively in (42). The exterior part was already analyzed by Collin &: Rothschild [16]. 
The subtracted power flow in (6) of the fields (42) has the radial dependence 

P T = ^ f E(r) x H*(r) • rr 2 dtt = (M) 
2 Jn 2|R^( K )|2 

in the exterior region r > a and vanishes in the interior region r < a. As the spherical 
vector waves are orthogonal over the unit sphere they can be analyzed separately. Their 
integrals are divided into its angular and radial parts. To simplify the notation, we in- 
troduce the normalized energies and w^j outside and inside the sphere, respectively. 
They are given by, see App. C for details 

«#= I \" ( iL(kr)\ 2 k 2 r 2 dn-ldkr 

J K JS1 



-~ " ^(1 h! 1 ^)! 2 - Re{h^(/c) h£»}) (55) 



for r = 1 and for the r = 2 modes 

"4? = I" ( |u^(fcr)|WdO-ldA;r = -Re{«h! 2) (K)(Khf 1) (K)) / } + ^ ) . (56) 

The corresponding normalized energy in the interior of the sphere is given by the 
integrals 



w® = f f \u { ^ ml (kr)\ 2 k 2 r 2 dndkr = [ x 2 \ j ; (x)| 2 dx 
Jo Jn Jo 

= y(j 2 («)-j i -i(«)j i+ i(«)) (57) 

and 

4i = f K [ |u^(fcr)| 2 fc 2 r 2 df]dfcr = -Re{Kj,( K )(Kj i ( K )/} + 4 ) . (58) 
Jo Jci 



We have the electric and magnetic Q factors 

qS»=|rWJ^ + -Jj^) (59) 

and 

<t, = |^W| 2 (^^ + ^|^). (60) 

VlRitWI 2 iRii'MIV 

After extensive simplifications we can rewrite them as 

oSuM-^ 51 ^ = «+«&.<«)• («u 
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Figure 2: Electric and magnetic Q factors for electrical surface currents J(r) = 
JoY T ami('f')5( r ~ a ) f° r ^ = 1,2. Power (solid curves) and far-field (dashed curves) stored 
energies. They differ by ka (61). a) TE (r = 1) modes, b) TM (r = 2) modes. 



where Q T -^ M denotes the electric and magnetic far- field type Q factors in (46), (47), (52), 
and (53). Note that the subscript EM is used to denote E and M in (61). The difference 
k = ka is consistent with the interpretation of a standing wave in the interior of the sphere, 
c/., (11). Moreover, the expressions (52) and (53) unifies the TE and TM cases and offer 
an alternative to the expressions in [20], here we also note a misprint in (6) in [20]. 

3.4 Numerical example for spherical shells 

The electric and magnetic Q-factors are depicted in Fig. 2 for I = 1,2. The relative 

differences are negligible for small ka where Q is large. For larger ka, where Q can be 

small, the relative difference is significant although the absolute difference is exactly ka. 

We also note that the Q factors oscillate and can be significant even for large ka. This 

is mainly due to small values of (ka) that can be interpreted as a negligible radiated 

power. Moreover, the Q-factors related to the far-field type stored energy (7) is negative 

in some frequency bands. The corresponding Q-factors related to (6) are always non- 

(P) (P) 

negative. Moreover, it is observed that Q n M > Q n E for low ka but has regions with 

(P) (P) 
Q\i,M < Qn,E for lar §er ka. 

To further analyze the negative values of (7), we depict the stored electric and magnetic 
energy density over spherical shells related to (7) in Fig. 3 for a TE spherical current sheet 
with radius ka = 1 and a coordinate system with origin at the center of the sphere. The 
results confirm that the far-field (7) and power flow (6) stored energy densities are identical 
outside the sphere. It is also seen that the far-field stored energy density is negative in parts 
of the interior region of the sphere, r < a, whereas the power flow stored energy density 
is non-negative. We also notice that the stored energy density is discontinues at r = a 
except for the far-field type stored electric energy. This is consistent with the boundary 
conditions that states that tangential components of the electric field are continuous. 
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normalized 




Figure 3: Illustration of the stored electric (E) and magnetic (M) energy densities for the 
TE (t = 1) mode generated by currents on a spherical shell with radius ka = 1. Power 
(solid curves-P) and far-field (dashed curves-F) stored energies. The energy densities are 
normalized with the radiated power and integrated over spherical shells to emphasize the 
radial dependence. The angular distribution is also depicted. 



4 Coordinate dependent term 



The stored electric (26) and magnetic energies contain the potentially coordinate de- 
pendent part W^S* defined in (28). Lets assume that = Wj F ^ f ™ ™ Q ™ 



EM,0 



for one co- 



ordinate system. Consider a shift of the coordinate system r — > d + r and use that 
r\ — r\ — > r\ — r\ + 2d ■ (r\ — r 2 ). This gives the coordinate dependent term 



w( F 2) _ w( F 2) , Uft . W 

EM,d - VV EM,0 + Ka VV ' 



(62) 



where W = W p + Wj and 
W 



2e JvJv 8vrA:ri2 




and we used (80), the identity 
sin(/cri2) 



p(ri)e 



-ikr-ri 



4vre 



dY 1 



ke 



,(r)\ 2 rdQ (63) 



Vi 



Airkr 



1.2 



ik lim / rG(\r - n\)G(\r - r 2 \Y dS 

r-+<x>J\ r \ =r 



ik 
16vr 2 



lim / r e -' ik Hri-r2) dQj ( 64 ) 



and the far-field potential (15). Similarly, the current part is 

W 3 = -^ [ [ J(r 1 )-J(r 2 rV 1 S |^dV 1 dV 2 = -- 1 -fe/ | A^rtfr dfi. (65) 
2 7y7y 8vrA;ri2 z 4/i 7 Q 



And with (18) totally 



W 



eo 
4/c 



|F(f)| 2 f dn. 



(66) 
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The corresponding Q factor is shifted as 

( f 2 ) _ -kd- J n r\F{f)\ 2 dn 

AQ ™ ~ "IpWdiT' (67) 

where we see that | AQ^ M \ < ka for all coordinate shifts within the smallest circumscribing 
sphere. 

Consider a spherical current sheet to illustrate the coordinate dependence. Let the far 
field be F ~ a\ YieOi +«2 Y20II, £-e., a combination of a z directed magnetic dipole and 
a y directed electric dipole. This gives ZiQj^ = —kx/4 and with a coordinate system 

centered in the sphere also Q^^o = as then r\ = r 2 giving Q^Md = —kx/4, where 
x = d ■ x and d is the vector to the center of the sphere. 



5 Small structures 

Evaluation of the stored energy for antenna Q is most interesting for small structures, 
where Q is large, e.g., Q > 10, and can be used to quantify the bandwidth of antennas [10, 
11, 15, 21, 22]. The low- frequency expansion of the stored energy are presented in [8-11]. 
Here, we base it on the low-frequency expansion J = j(°) + kJ {l) + 0(k 2 ) as k -»• 0, where 
V • = and the static terms J*- - 1 and po = — iV • /co have a constant phase. 

It is illustrative to compare the corresponding asymptotic expansions of the Q-factor 
components in (26). The coordinate dependent part vanishes if J and p(r) have constant 
phase. This gives 

Im{p(n)p*(r 2 )} = Im{(p (n) + fcpi(n))(p5(r 2 ) + kp\{r 2 ))} + 0{k 2 ) 

= klmipoi^pKri) + /9i(ri)p5(r 2 )} + 0{k 2 ) (68) 

as k — > and similarly for J. The different parts of the stored energy (26) contribute to 
the Q-factor asymptotically 

Q ™ ~ Qem ] - Ya> Q ™ ~ ka (69) 

as ka — > 0, where a is the radius of smallest circumscribing sphere and the coordinate 
system is centered inside the sphere. 

We can compare the expansion (69) with the Chu bound [15] 

gchu = (^F + ^' (70) 

where it is seen that Qchu has components that are of the same order as and 

and hence that these terms are essential to produce reliable results. This is also the 

conclusion from Sec. 3 in (61), where it is shown that the Q-factors differ by ka. 

The coordinate dependent part is negligible for small structures and of the 

same order as the difference between the far-field (7) and power (6) type as seen by 
the bound (11). We also note that the importance of Q diminishes as Q approaches unity. 
This also restricts the interest of the results to small antennas. The importance of 
for larger structures can however not be neglected. 
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6 Antenna examples 



6.1 Strip dipole 

Consider a center fed strip dipole with length I and width £/100 modeled as perfectly elec- 
tric conducting (PEC). The Q-factors (10) determined from the integral expressions 
in (27) and Q^p in (31), the circuit model [25], and differentiation of the impedance [21, 22] 
are compared in 4a. The circuit model is based on a the circuit representations of the 
lowest order spherical modes [26] with the lumped elements determined with the approach 
in [25]. The Q-factors from the circuit model approximates the integral expression very 
well for £ < 0.3A but starts to differ for shorter wavelengths where the circuit model is less 
accurate. The Q factors from the differentiated impedance is [21,22] 

where ' denotes differentiation with respect to uj and Z m is the impedance Z = R + jX, 
with j = — i, tuned to resonance with a lumped series (or analogous for lumped elements 
in parallel) inductor or capacitor 

Z m (u) = Z(u) - Jj*M^° if X M < ( 72 ) 
{jX^uo/u ifX(^ )>0. V ; 

In addition to the Q factor in (71), we determine the stored energy in the lumped element 
normalized with the radiated power as 

- W (73) 

giving the electric and magnetic Q factors 

Q m = \Q m ifx(c )<o 

^ E ifX(u, )>0 1 J 

and 

(Z)_/Q (Z) ifX(^ )>0 

ifXM<o, 

respectively. 

The difference — is also depicted in 4b. It is seen that the difference is 
negligible for the considered wavelengths. Curve (1) shows computed with a finite 
difference scheme. The curve is sensitive to noise and the used discretization. The noise is 
suppressed by approximating the impedance with a high order polynomial and performing 
analytic differentiation as seen by curve (2). 

6.2 Loop antenna 

The computed stored electric and magnetic energies for a loop antenna are depicted in 
Fig. 5. The loop antenna is rectangular with height £, width £/2, vanishing thickness, 
and is modeled as perfectly electric conducting (PEC). It is seen that the magnetic energy 
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Figure 4: Illustration of the Q factor for a center feed strip dipole with length £ and 
width £/100. The Q factors are determined from the stored energies (27) and (31) and 
from differentiation of the impedance (74) and (75). a) electric and magnetic Q-factors 
from (27), (31), the circuit model (dashed curves), and differentiation of the impedance 
Q( z ). b) difference between the computed Q-factors — Q^ z \ where is computed 
from a difference scheme and analytic differentiation of a high order rational approximation 
in 1 and 2, respectively. 



dominates for low frequencies. It changes to dominantly electric energy at approximately 
A « 6£ or equivalently A C/2, where C = 3£ denotes the circumference of the loop. 

In Fig. 5, it is seen that the Q-factors determined from the stored energies (27) and (31) 
and from differentiation of the impedance agree very well for Q > 10. The difference in- 
creases for lower Q values. This is consistent with the increasing difficulties to approximate 
the impedance with a single resonance model [22] and the potential ka ambiguity of the 
far-field stored energy (11). Here, it is also important to realize that the concept and 
usefulness of the Q-factor is increasingly questionable as Q decreases towards unity. 

6.3 Inverted L antenna 

An inverted L antenna on a finite ground plane is considered to illustrate the usefulness 
of the stored energies for terminal antennas. The antenna has total length i and width 
1/2, see Fig. 6. The electric and magnetic Q factors are depicted in Fig. 6. It is seen that 
Qgjjj and Q^M a S ree we ^ f° r Q ^ 10, that is for approximately £ > A/3 or below 1 GHz 
for 10cm chassis. The results start to differ for larger structures, where e.g., ~ 5 

and ~ 2 at £/X = 0.4 or ka ~ 1.4. For this levels of the underlying single 

resonance model [22] is problematic and hence reduce in accuracy. Moreover, 
as an approximation of Q have a relatively large uncertainty bound (11) for ka ~ 1.4. At 
the same time Q is low enough to be considered less useful as a quantity to estimate the 
bandwidth, e.g., Q 2 corresponds to a half-power bandwidth of 100%. 

7 Conclusions 

The analyzed expression (7) for the stored energy defined by subtraction of the far-field 
energy density from the energy density is mainly motivated by the formulation of Collin 
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0.5 1 1.5 



Figure 5: Illustration of the Q factor for a rectangular loop antenna with height I and 
width Ijl. The Q factors are determined from the stored energies (27) and (31) and from 
differentiation of the impedance (74) and (75). 



& Rothschild [16] and the expressions by Vandenbosch in [8]. We show that the stored 
energy (7) is identical to the energy in [8] for many currents. However, some current 
densities have an additional coordinate dependent term. This term is very small for 
small antennas but it can contribute for larger structures. Here, it is also important to 
realize that the classical definition [16] with the subtracted power flow (6) is inherently 
coordinate dependent. The identification of the energy expressions in [8] with (7) offers 
simple interpretation of the observed cases with a negative stored energy [11]. The analysis 
also suggests that the resulting Q factor has an uncertainty of the order ka. This is 
consistent with the use of the results for small (sub wavelength) antennas [11, 13], where 
ka is small and Q large. 

The energy expressions proposed by Vandenbosch in [8] are very well suited for op- 
timization formulations as they are simple quadratic forms of the current density. The 
quadratic form is very practical as it allows for various optimization formulations such as 
Lagrangian [11] and convex optimization [13] and has already led to many new antenna 
results. Their resemblance of the electric field integral equation (EFIE) makes the nu- 
merical implementation very simple. Analytic solutions for spherical structures show that 
the Q in [8] and [20] differ by ka and this is interpreted as the far-field in the interior 
of the sphere as seen from (7) and (6). The new formulation also produce simplified ex- 
pressions (61) the unifies the TE and TM cases for the Q factor of current densities on 
spherical shells [20]. Numerical results for dipole, loop, and inverted L antennas are also 
used to illustrate the accuracy of the energy expressions. 
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Figure 6: Illustration of the Q factor for an inverted L antenna. The circumscribing 
rectangle has height t and width £/100. The Q factors are determined from the stored 
energies (27) and (31) and from differentiation of the impedance (74) and (75). 
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A Green's function identities 

Multiply the Helmholtz Green's function for G\\ (V 2 + k 2 )G\ = —5(r — n) with G%, and 
similarly for G 2 . Adding the results together with a standard vector calculus identity gives 
2(VGi • VG* 2 - k 2 GiG* 2 ) = Gifa + G* 2 8 x + V 2 (GiG^). Integration yields the identity [8] 



/ 

Jm 



VG(|r - n|) ■ VG*(|r - r 2 \) - k 2 G(\r - ri\)G*(\r - r 2 \) dV 

co${k\r\ — r%\ 



(76) 



47r|ri — r%\ 

where we used Gauss's theorem together with the observation that 

V(GiG5) -»■ -fe ifcf " (r2 - ri) /(87r 2 r 3 ) (77) 

for large enough radius. 

The /c-derivative of the Helmholtz Green's equation for G\ is (S7 2 +k 2 )dkG\+2kG\ = 0. 
Similarly to the derivation of (76) we multiply with G 2 , and repeat the procedure with 
the ^-derivative of G 2 . Adding the result and applying vector calculus identities to move 
V 2 away from the ^-derivative results in the identity 

AkGiG* 2 = 5 2 d k Gi + 5id k G* 2 - V • q, (78) 

where 

r q = r- {G x Vd k G 2 - (d k G* 2 )VGi + G* 2 Vd k G x - (9 fc Gi)VG^) 



—k 
8vr 2 r 



\, ., „,1 



2 + - (f • (n + r 2 ) + i[|f x n| 2 - \r x r 2 \ 2 }) + 0{^)\ e -*Hn-r 2 ) (79) 
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for large enough radius. Collecting term of decay rate r _1 on the left-hand side and the 
remaining terms on the right-hand side. Integration over a large sphere, together with 
Gauss's theorem and elementary integrals results in 



/ 

Jm 



G(|r-ri|)G*(|r-r 2 |)- 
sin(/cri2. 



e -ifc(ri-T-2)-f 



167r 2 r 2 



dV 



8vr/c 



+ i 



8vrfc 2 7 



-(sin(/cri 2 ) — kryi cos(£xi 2 )) 



12 



sin(fcri 2 ) , . (n + r 2 ) • (n - r 2 ) . ,. . 

+ 1 ~ h{kri2) 



8tt/c 



87rri2 



sin(fcri 2 ) . (ri + r 2 ) sin(fcri 2 ) 
- 1 : • Vi- 



8tt/c 



8nkr 



12 



(80) 



Herej 1 (z) = (sin(z) — zcos(z))/z 2 and ri 2 = \r± — r 2 |. Note that (80) generalizes the result 
in [8] to the case T\ + r 2 ^ and shows that the integral depends of the coordinate system. 
It also shows that it is necessary to specify how the integration over R 3 is performed, i.e., 
here as the limit M 3 = {r : lim^^oo \r\ < ro}. 

B Spherical waves 

The radiated electromagnetic field is expanded in spherical vector waves or modes [24]: 



^Li(kr) = R^(kr)Y la mi(r) 



u 



laml 



(p), 



R^(kr) Y 2 ami(r) + #+1)^ Yami(r)r 



(p) 



(kr) 



U 



(p) 
3crml 



{kr) 



(p \kr) 



(kr) Yaml(r)r + ^T^T)*-^! Yw(^) 



where r is the spatial coordinate, r = r/r, r = \r\, k the wavenumber, and R 
the radial function of order I: 



(p) 



( M 



R#(«) 



= < 



AC 9« 



(81) 



(At) are 



(82) 



For regular waves (p = 1) zf^ = is a spherical Bessel function, irregular waves (p = 2) 

zf^ = n^ is a spherical Neumann function, and outgoing waves (p = 3) = h[^ is an 
outgoing spherical Hankel function. The indices are / = 1, . . ., m = 0, . . . , / , a = {e, o}, 
see [27,28]. In addition, Yramlif") denotes the vector spherical harmonics defined as 



Ylaml{r) = 



=V x (rY™ f (f)) 



(83) 



and Y2ami( r ) = r x Yiaml( r ) where Y a mi denotes the ordinary spherical harmonics [28]. 
There are a few alternative definitions of the spherical vector waves in the literature [24, 
27, 28]. Here, we follow [27] and use cosmcp and smmcp as basis functions in the azimuthal 
coordinate. This choice is motivated by the interpretation of the fields related to the first 
6 modes as the fields from different Hertzian dipoles. The modes labeled by r = 1 (odd 
v) are TE modes (or magnetic 2 z -poles) while those labeled by r = 2 (even v) correspond 
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to TM modes (or electric 2'-poles). We note that the derivatives of R™(k) are easily 
expressed in spherical Bessel and Hankel functions, i.e., 

dn ~ -Rg l(l + l)-K 2 (p) (M) 
I « + K 2 *, r-2. 

The Green functions are expanded in spherical waves to analyze spherical geometries. The 
scalar Green's function has the expansion [28] 

e ife|n-r 2 | 

G(n-r 2 ) = — , =ik^ ii(kr<)hl >(k r> ) Y am i(ri) Y CTmJ (r 2 ), (85) 

crmZ 

where r< = min{|ri|, |r*2 1} and r> = max{|ri|, |r2|}, and Y am i denotes the spherical 
harmonics. In addition, the full Green's dyadic, G = IG, can be expanded as [28] 

G(n - r 2 ) = ik £ "rL(kr<) u£„,(*r>), (86) 

rami 

where r = 1,2,3. We also use the frequency derivatives of the Green's function and the 
Green's dyadic expansions. 

C Volume integrals 

The volume integrals of the spherical vector waves are given by integrals of spherical 
Hankel functions as evaluated here. We have 

J x 2 z 2 p {x) dx = y (z 2 (x) - Zp-i^Zp+^x)) . (87) 

For the spherical Hankel function z p = hp 1 ' 1 = } p +i %, we have 

j x 2 \^\x)\ 2 dx = ^ (|h«0r)| 2 -Re{hp%(x)hpil(x)}) . (88) 
To evaluate the stored reactive energy outside a sphere, we need the result 

£° (x 2 \ h«(x)| 2 - 1) dx = a - | (| h«(a)| 2 - Rejh^a) h^(a)}) . (89) 
The corresponding internal energy is 

J\ 2 \^(x)\ 2 dx = | (j 2 (a) -j p _!(a)j p+1 (a)) . (90) 
We also have for r = 1 

/ |u lCTm Kfcr)| 2 dV= f \\^\kr)\ 2 r 2 dr (91) 
J[o,6]xn J a 

For r = 2, we use 

k\ U 2c rrnl | 2 = k U 2cr ml ■ u 2<7mZ = V X U lcrm 2 • 

= V • (u laml x u* CTmi ) + u laml -V x u* CTmi = V • (u laml x u* 2aml ) + fc| u lami | 2 (92) 
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and hence 



/ \u2aml(kr)\ 2 dY =1 

J[a,b]xfl fc 

where we have used the Wronskian relation z*z' — z'*z = —2i/x 2 , for z = h} 1 ^^), and the 
recursion relations for the spherical Hankel functions [4] in the last steps. The terms can 
be evaluated as b — > oo by considering the asymptotic behavior of the spherical Hankel 
functions. 



^ R (3) 

h (l)* K 2* r 2 

' kr 



+ I |h z (1) (Ax)|Vdr, 



(93) 
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